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A strong-coupling-perturbation theory around the Atomic 
Limit of the Anderson model with large U for a local- 
ized /-orbital coupled to a conduction-electron band is pre- 
sented. Although an auxiliary-particle representation is not 
used, application of the canonical Wick's theorem is pos- 
sible and yields an expansion in the hybridization V via 
dressed skeleton-Feynman diagrams. The Self-Consistent T- 
Approximation is constructed as a ^-derivable approxima- 
tion. From a numerical solution of self-consistency equations 
the /-electron-excitation spectrum is investigated. Compar- 
ison to the Non-Crossing Approximation is made in virtue 
of exact formal relations and numerical results. An exten- 
sion of this Feynman-diagram approach to the Anderson- 
lattice model is indicated, and application within the Local- 
Approximation scheme (limit of infinite spatial dimension) is 
given. 



I. INTRODUCTION 

Theories oLliighly correlated electron n systems like 
Kondo AlloysEla , Heavy-Fermion SystemsBu and High- 
T c Superconductors! 5 ! mostly start from model Hamil- 
tonians likcpthc Anderson or Hubbard model and its 
extensionstrEI . Recently a lot of interest has also beep 
taken in the multi-channel version of the Kondo modelE3 , 
to explain non-Fermi-liquid beha\aar.-abserved in non- 
magnetic 2-level-impurity systemstoxJ and certain U- 
based Heavy-Fermion compoundsolij . Common to 
these models is a single impurity or a lattice of impu- 
rities endowed with a local interaction. In Hubbard- 
like models considered here, electrons on localized d- 
or /-orbitals experience a Coulomb repulsion U . For 
the materials mentioned above U is the largest elec- 
tronic energy, and conventional weak-coupling theories 
appear not feasible. Besides several mainly numerical 
techniques, strong-coupling-perturbation theory around 
the limit of isolated localized orbitals (Atomic Limit) is 
appropriate for the problem with very large U . The mer- 
its of a self-consistent perturbational approach are con- 
tinuous excitation spectra calculated from Green's func- 
tions at the real frequency axis, and the possibility to 
separate quasi-particle formation from residual interac- 
tions. On the other hand, (diagram) rules are required 
which allow for a partial re-summation of the pertur- 
bation series to infinite order. In general this is diffi- 



cult to achieve in a strong-coupling approach, where the 
Coulomb-interaction operator ~ U is part of the unper- 
turbed Hamiltonian, whereas the small hybridization be- 
tween localized orbitals is taken as the perturbation. _ 

The earliest strong-coupling-perturbation approachll3 
deals with time-ordered Goldstone diagrams obtained 
from an expansion of the Hamiltonian's resolvenili3 with- 
out Wick's theorem and is well established for the 
Anderson-impurity model (for a Review see e.g. Refs. 
[ll||l9],|]) . There is nevertheless enduring interest in 
alternative approaches, since an extension to e.g. the 
Anderson-lattice iwdel imposes difficulties on the par- 
tial re-summationotHl , owing to the use of Goldstone 
diagrams, i.e. the loss of Wick's theorem. In Slave- 
Boson theories auxiliary (i.e. unphysical) particles are 
introduced together with a fluctuating field, which rep- 
resents the constraints . Conventional many-body tech- 
nique based on Wick's theorem can then be used in 
the calculation of fluctuation corrections to a Pjddie 
point, for impurity and also for lattice modelscjTEl. 
However, one has to start from the Mean-Field ground 
state (reviewed in e.gJtefs. |3l| , j32"| ), which exhibits bro- 
ken gauge symmetryo with an admixture of unphys- 
ical states which influence sum rulesOO . A variety 
of alternative diagrammatic approaches have been pro- 
posed, which start directly from the Atomic Limit rather 
than the Slave-Bose*i -Mean-Field ground state, like cu- 
mulant expansjjanaipEiJ , diagram-techniques for Hub- 
bard operatorsCjEj , and slave-particle methods where 
the Mean-Field phase is suppressed in keeping thq 
straint exact in each order perturbation theor) 
Other attempts-seek Wick's theorem in the limit of zero 
temperaturec3'EJ . 

The purpose of this paper is to present a self-consistent 
skeleton-diagram expansion around the Atomic Limit 
of the Anderson-impurity model with large U , which 
is built on Wick's theorem in its conventional formtj . 
Neither auxiliary degrees of freedom nor additional con- 
straints are introduced. The theory of lattice models is 
also considered in the light of-tbis new approach. 

In Sect, [n] an expansion! 5 in the hybridization 
V using conventional Feynman diagrams is re-derived 
and generalized. The Self-Consistent T- Approximation 
(SCTA) is constructed in Sect. Ill as a Conserving Ap- 
proximation for a set of fermionic Green's functions, and 
is specialized to the Kondo regime in Sect. |v|. Numer- 
ical results for the /-excitation spectrum are presented 
in Sect. [v|. Here also the relation to the Non-Crossing 
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Approximation is enlightened via numerical calculations, 
exact sum rules and spectral decompositions, and a cer- 
tain limit of the self-consistency (SCTA-) Equations. In 
Sect, fvj we turn to the Anderson-lattice model and inves- 
tigate the Feynman-diagram expansion on the lattice as 
well as the application of the so-called Local Approxima- 
tion. Re sults are summarized, and conclusions are drawn 
in Sect. |VIl|. Appendices [a| to |j5| contain calculational 



details not included in the main text. 

II. EFFECTIVE HAMILTONIAN AND 
APPLICATION OF WICK'S THEOREM 

A. Transformation of The Model Hamiltonian 

We start from the Anderson Hamiltonian! 

H = H C + H loc + H& + H[ 2 

of metallic conduction electrons H c coupled to a localized 
d- or /-shell H loc through a hybridization term + 

#12 : 



fe.cr 



(2.1a) 
(2.1b) 



= -±= E(l ~ ni a ) (Vkftcb, + h.c.) , (2.1c) 



nl = flf a , and Ngz denotes the number of wave vectors 
k in the 1. Brillouin zone. The one-particle level e£ of 
the localized orbital (designated as /-level and /-orbital, 
respectively) and conduction-electron band-structure e ka 
for spin a are measured relative to a fixed chemical po- 
tential. Conduction electrons form a symmetric band 
—D < £k < D with cut off D , and the /-level is situ- 
ated within the filled Fermi sea of conduction electrons, 
— D < el < . Indices k or a will be omitted in the 
following, if appropriate. 

The decomposition of the hybridization into and 
is a suitable starting point for strong coupling per- 
turbation theory, where the Atomic (Vk = 0) Limit 
H c + H loc is taken as the unperturbed Hamiltonian: The 
first part of the perturbation induces transitions only 
within the Fock space sector 7io,i which consists of all 
states with empty or singly occupied /-orbital. The sec- 
ond part H\ 2 mixes 7io,i and its complement TL2 contain- 
ing all states with double /-occupancy. The sector Ti 2 
involves the Coulomb repulsion U , and according to the 
hierarchy 



F< (|e / |,i?) < U 



it is convenient to eliminate at first H^ 2 in leading order 
V/U through a sequence of two canonical transforma- 
tions. The calculation is outlined in Appendix ^ and 
leads to the Hamiltonian 



H 



H" = H c 



■H- 



(2.3) 



In the effective interaction H J appearing here, an ex- 
change integral is simplified for large U and Vk = V to 
Ju :=-\V\ 2 /(ef + U)<0, 



N BZ tr^ 1 

k,q,a 



J u [flc\ _ a f- rT c k<J + flc[ _ ty c k ^ a f (T 



(2.4) 

The spin part of H J may also be written as — 2Jjjs c ■ 
and reveals the anti-ferromagnetic coupling of 
conduction-electron-spin density s c to the /-orbital spin 
s-^ due to virtual excitations into states with doubly oc- 
cupied /-orbital. 

In continuing the sequence of canonical transforma- 
tions we could also eliminate Hq[ and would reproduce 
the well known Kondo exchange Hamiltonian involving 
only spin degrees of freedom with a total coupling con- 
stant J tot = I V\ 2 /ef + Ju , equivalent to the result of Ref. 
|53| . We do not pursue this way, which is appropriate if 
the /-level lies well below the conduction band. 

An essential result from the above consi dera tions is 
the fact that the effective Hamiltonian Eq.(|]^) leaves 
Hilbcrt space sectors Ho,i and TL 2 decoupled; this enables 
the expansion in conventional Feynman diagrams given 
below. 



B. Wick's Theorem in Strong-Coupling Expansion 

For a study of the Anderson model's one-particle- 
excitation spectra and transport properties the quantity 
of interest is the /-Green's function 



F a {T-r') = -{T{U{r)ft{T')}) 



(2.5) 



The symbol T{. . .} orders canonical /-operators /£ ac- 
cording to imaginary time r, r' ; thermal expectation val- 
ues (. . .) and Heisenberg operators are formed with the 
Hamiltonian H" . 

By insertion of 1 = (1 — P 2 ) + P 2 , P 2 = n^n{ into 
Eq. (|2.5|) , the propagator is split into 



(2.2) 



F a (r - t') = F l r{r - t') + F^ h (r - r') , (2.6) 
with a low-energy part _F low and a high-energy part 

^high ^ 

F 1 ° w {t-t')= (2.7) 
-<T{ ([1 - nUfc) (r) ([1 - nL a ]fi) (/)}) , 

F^ igh (r-r')= (2.8) 
-(T{(n f _J a )(T)(n{ a rt)(T>)}) . 
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The projector P2 onto subspace H.2 commutes with H" , 

[P 2 ,H"] = 0, (2.9) 

i.e. double /-occupancy is con served, and therefore mixed 
terms not written in Eq.( |2.6|) va nish. 

The low-energy part Eq. (|2.7|) contributes to the one- 
particle-excitation spectrum fr(u>) of /-electrons in the 
frequency range < < D we are interested in, 
whereas the spectrum of F hlsh is well separated from 
that by U , due to excitations into H2 ■ Thus the high- 
energy part will be ignored in the following and we take 
F = F low , that is 



F a (r-r') 



(2.10) 



_2_Tr[e-^"T{ (/- CT /! CT / CT )(r) (f-afUt) (/) }] 

with the partition function Z = Tr[e~P H } and the in- 
verse temperature [3 — 1/ksT . 

Our goal is a perturbation expansion of the /-Green's 
function in the hybridization and the exchange cou- 
pling term H J appearing in Eq.(2.3). In order to ap- 
ply Wick's theorem on the course, two points have to be 
noted: First, on account of the conservation law Eq.(|2.9|) 
the Fock-space sector H2 involving double /-occupancy 
is completely decoupled from the dynamics introduced 
by the perturbation. Second, states from TL2 do not con- 
tribute at all to the trace in Eq. ( 2.10D , owing to Eq. ( |2.9| ) 

and P 2 (f-afl a /P) = (/- CT /1 < 3 TJ )P 2 = 0. Thus the 
interaction term in the local /-Hamiltonian H loc is com- 
pletely ineffective, and we may assign any value to the 
parameter U in H loc , or ignore this interaction term for 
convenience: H loc — > H loc = s^n^ . Accordingly we 
may work with a modified Hamiltonian H" — > H , 



H = H {0) + H, 



01 



H J 



(2.11) 



where the unperturbed part is now bilinear in band- 
and /-operators, 

= e k J ka c ka + £ e f jlh ■ (2.12) 

Introducing properly normalized expectation values 

(...) = ^Trie-P" . . .] , Z = Tc[e- f>S ] , (2.13) 

z 



the one-particle Green's function Eq.( 2.10|) now reads 



Fa(r-r') = -F a (r-r'), 



(2.14) 



with 



F CT (r-r')= (2.15) 

(f-afUa) W (/-<x/l ff /l)(r')}T ■ 



In F the thermal average as well as time dependencies of 
(composite) operators A(t) = e TH Ae~ rH are determined 
by the modified Hamiltonian H . Note that double /- 
occupancy does contribute to Z and therefore Z differs 
from the partition function Z . 

An expansion in connected Feynman diagrams by ap- 
plication of Wick's theorem and linked-cluster theorem 
is now apparent, since F is a 3-particle Green's func- 
tion of a fermion system, with an unperturbed Hamilto- 
nian H and a tw o-part icle interaction + H J given 



through Eqs. Q2.lcD and fl2.4| ) . Following e.g. Ref. |54] a 
diagram of given order in V and Ju is constructed from 
bare fermionic Green's functions and the vertices shown 
in Fig. 0: A dashed line represents a local fermionic 
Green's function 



ILUl) 



1 



(2.16) 



with Matsubara frequency loi — (21 + 1)tt/ (3 . A full line 
corresponds to a band electron at the impurity site, with 
the fc-dependence of included, 



1 



N B z 



E 



\Vk\ s 



V 2 iuJl - Ska 



(2.17) 



For counting the order in perturbation theory a real av- 
erage V has been introduced, 



V 2 := 



1 



N B z 



(2.18) 



Vertices originating from the correlated hybridization 
are shown in Fig. |l| (a) and contribute a factor V 
each, those corresponding to H J carry a factor Jjj , and 
are depicted in Fig. [l] (b) . All vertices conserve direction 
of arrows (i.e. particle number), spin a and Matsubara 
frequency u>i . Each internal free frequency to be summed 
up is accompanied by a 1/(3 = ksT and an unusual con- 
vergence factor exp(— iujiQ + ) : its exponent differs in sign 
from standard rule s, ow ing to non-normal ordering of /- 
operatorsE3 in Eq.(2.1c) . Sums over internal wave vec- 
tors k are always independent and already included in the 
local band-electron line G c . The sign (— l") n °+Xext+nj f 
a diagram is determined as usual by the number n c of 
closed fermion loops, the ordering Xext of external time 
labels, and the number nj of vertices ~ Jjj . The sign 
contribution from vertices ~ V is compensated by the 
non- normal ordering (1 — w) = ff* of /-operators in 
Eq.Q. 

An interpretation of vertices and diagrams in physi- 
cal processes as indicated in Fig. [l] is obvious, but has 
to be considered merely as illustrative: The dashed line 
Qf(o) does not represent a true bare /- particle. 111 
the unperturbed limit V = , which is given byu 

F^^ l ) = (l-nl a )^G^^ l ) ; 
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hence the propagator G^ ^ lacks the reduced spectral 
weight factor of . Although this is only a minor 
effect here, differences become much more pronounced if 
full propagators are considered. 

For the remainder of this paper we turn to the U — * oo 
limit, that is, vertices ~ J\j in Fig. [l] (b) will be dropped 
according to the Hamiltonian H =>■ + . The 

partition function Z is now given through 



U -> oo : Z = Tr[e 



-/3H" 



] -> Tr[e^ K (l - ^ 



Z-Z c exp(-/3V4) . (2-19) 



Thus Z is accessible via an expansion of ln(Z) in linked 
diagrams. Z c denotes the partition function of bare con- 
duction electrons. 

The Hamiltonian Eq. ^jj) 
i.e. J = is well knowr 
ied using expansions in the hybridization via time 
ordered|-(jCk>ldstone) diagrams based on the Resolvent 
Method or equijsiajgit schemes involving constraints 



m the limit U — > oo , 
It has widely been stud- 



on auxiliary particlesoO . Nevertheless the possibility to 
apply the canonical Wick's theorem in a strong-coupling 
expansion has not been realized so far. 

Furthermore it is interesting to note that the approach 
proposed here is easily transferred to the spin-1/2 Kondo 
modeEl and multi-channel Kondo mode£3 , giving rise 
to a Feynman-diagram representation of the conduction 
electron's T-matrix: The additional states introduced in 
writing the local spin in fermions do not contribute to 
the T-matrix; thus any 'slave-fermion' constraint is un- 
necessary for a spin-1/2 impurity model, as has already 
been argued in Ref. [38|. 



III. SELF-CONSISTENT APPROXIMATION: 
GENERAL ASPECTS 

In the diagrammatic expansion of the 3-particle prop- 
agator F , which co vers t he dy namic s of one-particle /- 
excitations via Eqs.(2.14) and (2. IS), we encounter the 
logarithmic divergences substantial for formation of the 
Kondo effedH . This can already be confirmed for the 
simplest ladder diagrams shown in Fig. EjL which con- 
tribute to F . Details concerning representations of F will 
be given later on. Consider e.g. the particle-particle (pp) 
diagram "f pp {ivk) ■ Its analytic contribution at bosonic 
frequency iv\ = 2kir/(3 is 



7»»(u*) « V 2 p c (0) [ D de . /(£) , 
J-D lv k - £ - £ T 

where ksT <§C \e*\ and f(e + iuk) = /(e) denotes the 
Fermi function. Spin degeneracy has been assumed, and 
a flat conduction-band density p c {e) = — ^lmG c (e + 
i0 + ) = jp<d(D — |e|) inserted. After continuation to the 



real axis, iuk — > w±iO+ , the real part Re7 pp (w±i0+) de- 
velops for T^Oa logarithmic singularity at the /-level 
u> = e> , Similar results are obtained for the particle-hole 
(ph,ph) diagrams in Fig. 0. 

The strategy to be pursued is quite common to 
known perturbational strong-coupling approaches^ . Di- 
agrams showing logarithmic divergences have to be re- 
summed, whereby a small Kondo-energy scale .in the 
form T A = Dexp(e f /2V 2 p c (0)) is establishedBE] . Fur- 
thermore a renormaliaatiou—of propagators in skeleton 
diagrams is necessaryETE3'E£l , which leads to a smooth 
Abrikosov-Suhl resonance near the Fermi edge in the 
/-excitation spe&tsum and turns Ta towards the exact 
Kondo energyOL3 Tk ■ Accordingly we aim at an ap- 
proximation for the propagator F via summation ,-af a 
certain class of dressed skeleton-Feynman diagramsEJ to 
infinite order in V , involving the diagram elements shown 
in Fig. §. 



A. Skeleton-Diagram Expansion 

To begin with, consider the construction of proper self 
energies. Since the vertices from Fig. |l| (a) present in 
the U — > oo problem allow for connected diagrams with 
external lines of arbitrary type, i.e. G c (c-line, given in 
Eq.(pT7|)) or &^ (/-line, see Eq.(§Ig)), we get four 
self energies designated as 



(3.1) 



The type of outgoing and incoming line is indicated in 
the superscript. As usual, external lines are removed and 
irreducible (proper) self-energy diagrams cannot be split 
by cutting only one line of any type. Therefore not only 
a renormalization G^ * 1 — > & and G c — > G c takes place, 
but also mixed Green's functions G^ c and G c f emerge. 
The G aa are depicted as the various kinds of double lines 
shown in Fig. ^| (a) . In the imaginary-time domain, these 
propagators are defined as 



Gl(r-r')= -(r{/ CT (r)/t(r')}T, 

Gj c {r-r')= - (T{f a (T)Ci(T')}) , 

G c J{r-r')= - (r{a(r)/t(/)}[, 

G c Ar-r')= - {T{C a {r)Ci{r')}) , 

with local operators and = -?A — 



(3.2a) 
(3.2b) 
(3.2c) 
(3. 2d) 

i r) ,(t) 

k V ka ■ 

Green's functions Eq.( |3.2j ) and corresponding self ener- 
gies Eq^SJ) are now collected in a matrix propagator 



Gcr(iuJi) 

and matrix-self energy 
£ CT (wj/) = 



qi Gp 

G C J G°„ 



(3.3) 



(3.4) 
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The bare (V — 0) matrix G*- -* contains only diagonal el- 
ements G-fW and G c since mix-Green's functions G^ c ' c ^ 
vanish in the Atomic Limit. Full and bare matrix-Green's 
functions are shown in Fig. || (b) , followed by Dyson's 
equation in Fig. |^ (c) . Diagram rules for matrix propa- 
gator and vertex function as well as the vertex in matrix 
formulation follow from re-writing the perturbation 
in 'spinors' ty a = (f a , C„) , f CT = (/] , G+) . The result- 
ing vertex is displayed in Fig. || (a) . At last, we combine 
the vertex' two possible orientations as shown in Fig. 
^ (b) . The net result is, that diagram rules for unphys- 
ical 'free energy' ln(Z) (see Eq.( |2.13| )), Green's function 
G , self energy S and all vertex functions are quite anal- 
ogouS|-tp the interacting electron gasE.3 , with an unusual 
factorEa exp(— iuj n 0+) for each internal frequency ui n . 
Dyson's equation shown in Fig. ^ (c) is written 

(G^)r 1 = (G^or 1 - 

and is solved by inversion of a 2 x 2-matrix. The elements 
of Go- obtained this way are noted for later use, 



61 



i 



[G«°>] 



with 



Gl c = Gl 



G C J = 



H = W 



(3.5a) 



\G- 



K c [G 
-Gl = Gl 



/(0)i 



/(O)l-l 



Gl 



1 



\G%\ 
with 



l _ S c 
= 



[G; 



£c/ E /c 



G£ (3.5b) 
(3.5c) 
(3.5d) 



[G 



/(0)i 



The 3-particle propagator F a of Eq. ( [2.15] ) , which is 
the quantity that furnishes the /-excitation spectrum, is 
conveniently expressed thro ugh the self energies appear- 
ing in Dyson's equation Eq. l \>.rk : In close analogy to the 
well known T-matrix relatiorEa exact equations of mo- 
tion are utilized to show that 



G%(iu>i) = G c a (m) + G£(»Wi) V 2 F (7 (iuji) G% (iw,) , (3.6) 
with bare and renormalized conduction electron Green's 



function G% and G C G , as defined in Eqs. (|2.17| ) and ( |3.2d| ) 
respectively. In combination with Eq.(3.5d) we get the 
exact expression 



V 2 F a (iuo l ) 



l-G-im^iiui) 



(3.7) 



Alternatively this result is derived from considering the 
diagram series of F a directly. Apparently V 2 F equals 



the improper self energy of G c . Since F stands for a 3- 
particle propagator, it contains diagrams reducible with 
respect to band lines G c (i.e. which can be split by cut- 
ting a single G c ) compliant to particle conservation, in 
contrary to any two-particle Green's function. 



B. Mean-Field Theory and Infinite Ladder 
Summation 

The formalism developed so far serves as the basis 
for defining self-consistent approximations: An expan- 
sion of the matrix-self energy X in skeleton diagrams 
with matrix-type lines and vertices leads to consistent 
expressions for its elements E a a , a' a = ff,cf,fc,cc. 
These fully determine the self-consistency problem via 
Dyson's equation Eq.(|3.5|) as well a s the one-particle /- 
Green's function via Eqs. (3/7) and (£33)-. As a starting 
point we take the exact representationOLj for S given in 
Fig. ^ . To prevent over-counting the lowest order contri- 
bution ('Hartree diagram') is written explicitly, whereas 
all highetjorders are absorbed into the irreducible vertex 
functionO T 2 P . 

At mean-field level only the 'Hartree diagram' in Fig. 
H is taken into account, which leads to a self energy £ 
independent of frequency ui . The self -con sistency equa- 
tions resulting from Dyson's equation ( |3.5| ) can be solved 
analyticallyO at T — > . They possess, however, only a 
trivial solution: 



S^(iw,) = 2V(/! ff G_ ff ) 
5#(tw,) = ££ c ^) = V(l 



0, 



o, 



That is, the unperturbed (V = 0) Atomic Limit is re- 
covered in mean-field theory. This result fits into the 
discussion of the Kondo regime given in Sect. |v| below: 
The lowest order diagram (Fig. |5| left) considered here 
belongs to a whole class of skeleton diagrams which are 
negligible in the Kondo regime. 

In going to higher orders, attention has to be payed to 
symmetry relations fulfilled by the vertex function T 2 P : 
In Fig. [| the self energy S has been given in a way that 
T 2 p enters X in the particle-particle (pp) representation 
Tp(l', 2'; 1, 2) , with arbitrary labels attached to incom- 
ing (1, 2) and outgoing (1', 2') lines. Correspondingly 
it can be analyzed in the (pp)-channel, i.e. re-written 
through the Bethe-Salpeter equation 

rf,(l',2';l,2) =r»(l',2 / ;l,2) (3.8) 
+ r^(l',2';T,2)G(T ) T')G(2,2')r|(T' ) 2';l,2) , 

involving a kernel T pp which is two-particle irreducible 
only in the (pp)-channel. An integration J dl d2 dl d2 
is implicit, where e.g. J dl is shorthand for L dr\ 
— ±1) — /> c )- However, the self energy 

could also be written with reversed arrows on the spin 
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(— a) loop in Fig. 0, so that Tp would be replaced by 
its particle-hole (ph) representation r 2 (l', 1; 2, 2') . Now 
two alternative Bethe-Salpeter equations are apparent, 
where V 2 is analyzed in the (ph)- or (p/i)-channel with 

respective kernels T ph and T ph , 

r 2 (l', 1; 2, 2') = r**(l', 1; 2, 2') (3.9a) 
+ r^(l', 1; 2, 2') G(2, l')G(I, 2') r 2 (T',T; 2, 2') , 



r 2 (lM;2,2')=r"(lM;2,2') 



(3.9b) 



+ r ph (l',T; 2, 2') G(T,T)G(2, 2') r 2 (T', 1; 2, 2') . 

Additionally an interchange of the two incoming lines 
of rp in Fig. |5| or, if the particle-hole representation has 
been chosen, of T 2 does not alter £ , in virtue of the corre- 
sponding anti-symmetry of Tp , T 2 . In combination with 
the strict equality of particle-particle and particle-hole 
representations this anti-symmqtey is expressed through 
the so-called Crossing RelationsEl 



-rp(l',2';2, 



1) 



(3.10) 



= T 2 (l', 1; 2, 2') = -r 2 (l', 2; 1, 2') 
These should be fulfilled by any approximation. 



A ladder-type approximation to T 2 p ^ may be obtained 



Y ph 



from Eqs. (fj|) and with the kernels T pp , T ph 
each replaced by the bare vertex shown in Fig. ^ (b) . 
An iteration of the Bethe-Salpeter equations then yields 
three different ladder approximations for the vertex func- 
tion, with diagrams from only one of the three inequiv- 
alent channels included. Taking solely one of them into 
account, with reference e.g. to RPA, does certainly not 
comply to Crossing Relations. It is more suitable to add 
all three ladder sums and to remove the lowest order con- 
tribution (i.e. the bare vertex) twice afterwards, since it 
has been counted three times. The resulti ng ve rtex func- 
tion rp = r 2 fulfills Crossing Relations ( 3.10 ) and con- 
tains no over-counted diagram. Furthermore it includes 
all types of loops from Fig. ||to infinite order. Eventually, 
insertion into Fig. |^ leads to the skeleton approximation 
for the self energy £ displayed in Fig. [|. 

The ladder summation obtained is considerably sim- 
ple, since solely bosonic propagators representing ladder 
sums have to be treated self-consistently. Ar«pije-system- 
atic approach based on Parquet equationscSEjO would 
involve 'true' retarded vertex functions. Additionally the 
approximation is ^-derivable, i.e. a functional $ can be 
foundO with the propertyO 



(S CT (r'-r)) c 



(3.11) 



Thus all self energies S a a , a' a — //, cf, fc, cc can be 
obtained by variation of one functional $. As a conse- 
quence the ladder approximation is thermodynamically 

consistcntL-3 , i.e. expectation values (...) with a tilde are 



safely computed from corresponding Green's functions. 
Since the diagrammatic representation of $ xesembles 
that of the T-matrix theory of the electron gasEZl the self 
energy Fig. ^ will be referred to as the Self-Consistent 
T- Approximation (SCTA). 



IV. SELF-CONSISTENT SCHEME AND 
EXCITATION SPECTRUM IN THE KONDO 
REGIME 

So far an approximation for the Anderson model in 
the strongly correlated limit U — > oo has been derived, 
which does not rely on any specific parameter range such 
as the Kondo or mixed-valence regime. In the following 
the Kondo regime will be considered further, where the 
coupled self-consistency equations giv en through the ma- 
trix X and Dyson's equation Eq.( |3.5[ ) may be simplified 
using a non-perturbative small parameter. 



A. Small Parameter 

The small parameter available in the Kondo regime 
is given by the c-number Z/Z , which relates physi- 
cal quantities like the one-particle /-Green's function F 
to i ts co unterpart indexed by a tilde, as introduced in 
Eg. (PI: 



F a {T-T') = -F a (T-T> 



(4.1) 



At U — ► oo the prefactor is given through Eq.(2.19) , 

Z/Z = 1 + (Z c /Z) cxp(-2/3e / ) . (4.2) 

We estimate the order of magnitude of the partition func- 
tion Z via 

Z/Z c = exp(-j9(n - il c )) » exp(-(3(E G - E c )) , 

where the free energy Q has been replaced at low temper- 
ature T — > by the ground-state energy Eq . Quantities 
Z C ,Q C , E c correspond to bare conduction electrons. As 
is well known, hybridization lowers the Atomic Limit's 



(V = 0) ground-state energy Eq 



(0) 



E € only by 



the small amount of the Kondo energy Tk to Eq w 
E { q ] -T K . Thus the estimate Z/Z c w exp(-(3(s f -T K )) 
holds, leading to 



- = l + i«i )S: = exp(/J(e' + T K )) 

Z K K 



(4.3) 



This provides a non-perturbative parameter k <C 1 , 
which is exponentially small in the Kondo regime < , 



T K « \e f \ , k B T « \e f \ . From Eq.(0) it now follows 



that F is of very small magnitude ~ (k) 1 (or smaller), 



G 



with F being an intensive thermal expectation value 
~ (k) = 1 (or smaller) . 

In a similar way information is gained on the mix- 
Green's f uncti ons G c} ' ^ c : Following the arguments given 
after Eq.( 2.10| ) , the relation 



(T{a(r)4(r')}) = |<T{C a (r)/t(r')}) 

holds, together with its complex conjugate involving f a 
and C\ defined belo w E q.(3.2) . Thus the off-diagonal 
elements of G in Eq.(|3.3[) are of order ~ (k) 1 , 



(4.4) 



GI c {t~t')~{k)\G c J{t-t')~{kY 



It is also useful to consider directly the spectral weight 
of F a , (1 — n£) = (Z/Z)(l - n£). It follows that 

dLjpf(u)[l-f(u)]~ K) (4.5) 

- OO 

with the Fermi function f(tu) and the spectral function 



#(a;) = -ilmG£(w + iO + )>0 



(4.6) 



Since the integrand in Eq.(4.5) is positive, the function 

^MtWMl-OO 1 (4.7) 

may serve as a small parameter, too. This property can 
also be derived from the fact that the spectrum p' de- 
velops a threshold behavior at energy E Th + Tk 
when the temperature goes to zero: In Appendix |D| it 
is shown that p f (u) ~ @(E Th - u) at T = . It is ex- 
pected that pf (u>) shows a singularity as lu approaches 
the threshold E Th from below. Nevertheless, at finite 
temperatur e th e singularity is removed, and the left hand 
side of Eq.(4.7) shows its maximum at uj « E Th with a 
value PiE^fexp^E 711 ) w k < 1 . 

With the small quantities Eqs.(4.4) and ( }4.7[ ) at hand, 
self-energy diagrams of leading order in k can be sepa- 
rated from those of higher order in k , which will be ne- 
glected in the Kondo regime. For that purpose we return 
to the explicit notation using the vertices and Green's 
functions shown in Figs. [I] (a) , | (a) , and dissolve each 
matrix-type diagram from Fig. 6 into its contributions to 
E a a , a' a = //, fc, cf, cc . Wher eas mix-Green's func- 
tion s co ntribute directly via Eq.(4.4) a factor k each, 
Eq.(4.7) is used within a loop theorem proven in Ap- 
pendix^]: A closed path in a skeleton diagram to E a a , 
which runs exclusively on equally directed Green's func- 
tions (double dashed lines), is sufficient for a prefactor 
k to that diagram. Vertices may be passed in arbitrary 
fashion from any incoming dashed line to any outgoing 
dashed line. Accordingly a closed path in general involves 
lines of either spin direction. As an example two dia- 
grams are shown in Fig. ^| (a) on the top and bottom, 
which possess one and two non-overlapping closed paths 
respectively. By application of the loop theorem these are 



of order ~ (k) 1 and ~ (k) 2 . Closed paths are indicated 
by dotted lines; in the diagram on the top only one of 
two equivalent paths is marked as an arbitrary choice. 

The resulting subset of ladder diagrams to be con- 
sidered in the Kondo regime is displayed in Fig. for 
the diagonal elements of E , i.e. E-^ and E cc . f is 
~ (k)° = 1 because no closed path as described above 
can be found in the diagrams shown in Fig. |§| (top), nor 
are any mix-Green's functions G^ c or G c ' present. The 
diagrams of (pp)-type not shown here as well as all (ph)- 
and (p/i)-diagrams to E" are at least ~ (k) 1 by applica- 
tion of the loop theorem , i.e. vanishing small in the Kondo 
regime. Fig. |7] (a) displays examples. Also any diagram 
containing at least one mix- Gree n's function G? c or G°* 
is ~ (k) 1 or smaller by Eq.( |4.4D and can be neglected. 
Thereby all diagrams of odd order V 2n+1 are unimpor- 
tant, as illustrated in Fig. [?] (b) . E cc on the other hand 
is in leading order ~ (k) 1 : Every diagram of the (pp)- 
type ladder sum indicated in Fig. || (bottom) shows one 
closed path on equally directed dashed lines. All (ph)- or 
(p/i)-type ladder graphs to E cc not shown here contain 
at least two non-overlapping closed paths. These con- 
tribute at order ~ (k) 2 and are omitted, as well as all 
diagrams containing at least two mix-Green's functions 
Qcf/fc ^ F ur thermore, it can be seen that all remaining 
ladder graphs to E cc with only one mix-Green's function 
contain in addition at least one closed path on dashed 
lines and thus are also negligible. 

The above-mentioned simplification of the self- 
consistency problem posed through Dyson's equations 
Eq.( |3.5|) follows from the fact that E-^ is ~ 1 in lead- 
ing order k , whereas E cc and /-c-mixing self energies 
are small, 



E^ ~ 1 , E c/ - E /c - (k) 1 , E cc - (k) 1 



(4.8) 



This is shown in A ppen dix |b] to hold in general for the 
Kondo regime. Eq.(3.5) now reduces to 



(4.9) 



and G% = G% - 1 , G ! a c - G°J ~ (k) 1 . Therefore & is 
renormalized self-consistently through E" given in Fig. 
^, whereas the conduction-electron Green's function G c 
remains bare. The magnitude of mix-Green's functions 
Qfc/c} i s reproduced consistently. The latter need not 
be considered in the Kondo regime. 

The local /-Green's fun ction follows from Eqs.(f.l) 



and (|3.7|) , (|3.5cQ with (^8|) as 



Z V 2 



(4.10) 



Although E cc ~ (k) 1 is unimportant for the renormaliza- 
tion of unphysical propagators, it determines the physi- 
cal /-Green's func tion F , which in combination with the 
prefactor Eq.(4.3) is of order 1. The property of E being 
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<l>-deriva ble st ated at the end of the last section, together 
with Eq.jll0|) takes the form 



d<5>\GfG c 



dGl{iuji 

This corresponds to analogous < 
Resolvent-Perturbation Theor\ 



F„{i<jJi) 



Z <9$[G-f,G c ] 



ZV 2 dG-(iu l ) ' 
ressions derived within 



B. Self-Consistency Equations 

In order to derive analytical expression for self ener- 
gies, two 'ladder elements' are defined as shown in Fig. 

7?(«Vfe) can be understood as an effective bosonic prop- 
agator with external frequency v k = 2kir/[3, which is 
further renormalized via a 'self energy' a(ivk) to 



ri(wfe) 



ir(iv k ) 



1 - n(iuk) cr(iv k ) 



(4.11) 



Spin degeneracy is assumed here and in the following, 
although II is spin-independent in general. The formula 
corresponding to Fig. (top) then reads 



lUJl) 



(4.12) 



with external and free fermionic frequencies io\ and Lu n — 
(2n + l)7r//3 respectively. An overall sign results from 
the single closed fermion loop in every diagram. With 
t he lo cal DOS of conduction electrons (see Eqs.(2.17) and 



\V,\ 2 



Nbz^ V 2 



5(e - e k ) 



(4.13) 



and a spectral function of the effective propagator II , 

p n (e) =-ilmn( £ + zO+) , (4.14) 



the frequency sum in Eq.(4.12) is performed through a 
contour integration, leading to 



(4.15) 



dede' p c {e)p n {e') 



[l-f(e)] + [-l-g(e')} 



iui 



■e — e 



A Bose function g(e') appears through f(e' — iu>i) = 
—g(e') ■ In the numerator a 1 — 1 = has bee n added, in 
order to utilize a property similar to Eq.(h7) , 



^( £ )[-l- 5 ( E )]~(K)i 



(4.16) 



Thereby the second term in Eq.(4.15) is of order ~ (k) 1 
and will be neglecte d in f avor of the first term ~ 1 . 

The estimate Eq.( 4.16 ) follows from considering the 
bosonic Green's function 



CT{(/T/i)(r) (/iVtV)}) 



(4.17) 



It acquires the form Eq.(hll) when its diagrammatic is 
considered in SCTA for the Kondo regi me. In a fash- 
ion similar to the argument following Eq.(h4) we relate 
this propagator to a suitably chosen thermal expectation 
value, namely 

((l-nf)(l-n{)) = (Z/Z)fi((r-^) = 0+) . 



With Eq.fl4.3j) this results in 

J dejP(e)[-l-g(e)]=ti(0 + )~(K) 1 



positive 



and Eq_dhl6) follows since the integrand is always 
In Appendix ^ this conclusion is drawn di- 



rectly from Eq.( [t.li ) . Similar to pf the spectrum p^{uf) 
shows a threshold behavior at uj = E Th z* + , 
as is shown in A ppen dix At finite temperature this 
also implies Eq.(4.16), following the arguments below 

Eq -& 

T he se lf energy is given through the first term of 
Eq.( 4.15 ) which reads after continuation to complex half 
planes iuji — > z , Im(z) 7^ , 

Y,ff(z) = V 2 { dep c (e)[l- f(e)]U(z + e) . (4.18a) 



This has to be inserted into G? given through Eqs.(^ 
and (frf), 



G f (z) =[z- e f - T, ff (z)] 1 . (4.18b) 
The right hand side of Eq. ( [4. 18a| ) requires 

U(z)=[n(z)- 1 -cr{z)]- 1 , (4.18c) 



according to the definition Eq.( 4.11 ) . it and a given in 
Fig. H complete the set of equations, 



a{z) = N,jV 2 / dep c {e)f(s)G f (z-e) 



(4.18d) 

with spin-degeneracy Nj = 2 , and 

%{z)= I dep f {s)[2f(e)-l]G f (z-£) . (4.18e) 



The latter two exp ressions follow in leading order ~ 
(k)° from Eqs.(|Cj) derived in Appendix |Cj. The 
self-consistency equations ( 4.18Q represent the Self- 
Consistent T- Approximation to the spin-degenerate U — > 
00 Anderson model in the Kondo regime. Parameters are 
the /-level e ' , hybridization matrix-element V , the con- 
duction band's density of states p c (e) , and the tempera- 
ture via the Fermi function /(e) . The Green's function 
& , or its spectrum pf (e) is by construction the inde- 
pendent variable to b e determined through a (numerical) 
solution of Eqs.fl4.18|) . 
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C. Physical /-Spectrum, Sum Rules 

For calculating the /-excitation spectrum we now turn 
to S cc : From inspection of Fig. || it follows 

J J ' ' tLui+e-e 1 



quite analogously to Eq.( 4.15 ) . Here, with Eqs.([4~7|) and 
(4.16) both terms of the right hand side are ~ (kt, and 
the resul t fulfills the general property Eq. ( |4.8| ) . With 
Eq.(4.10) the /-Green's function follows immediately, 
and the /-excitation spectrum reads 



(S(u) = - ilmF(cj + iO+) 



(4.19) 



de 



=p f {e)r^ + e) + pi{e)T^ + e) 



Projected Spectra have been introduced here, which are 
~ (k) and positive for all energies, 



f(s) = -pf(s)[l-f(e)}, 
f 1 (e) = ^p n (e)[-l-g(e)\ 



(4.20a) 
(4.20b) 



—f 

It is suitable to use an extra set of equations for p and 
p n , in analogy to the Defect Equations in NCA-theoryEl . 
These are derived in Appendix ^| and read 

f(cj) = G f (u) 2 V 2 f dep c (e)f{e)f(o J + s) , (4.21a) 



f(u) = fi( w ) NjV 2 I de P c {e)[\-f{e)]=p f {u-e), 

(4.21b) 



with fixed spin degeneracy Nj = 2 . 

For the spectra occurring in Eqs.( 4.1g| ) and ( 4.19 ), 
. 2 l h a set of exact sum rules is proven in Appendix 
Here we quote the results: The integrated spectral 
weight of JF and ff 1 is given by 



J dep f (e) = 1 , 

def l {e) = l-2 I dspf(e)[l-f(e)} 



(4.22a) 
(4.22b) 



with Eq.(4.5) the weight of pP can be set equal to one. 
Projected Spectra are related via 



de 



2f(e)-f(s) 



(4.22c) 



and the physical /-spectrum is normalized according to 



J dep f (e) = (1 - n f u ) = J dep f {e) . (4.22d) 

Note that the integrated /-electron spectral weight is not 
equal to unity, since excitations at high energies ~ U are 



excluded here; compare to Section [I B . 

In Appendix H the SCTA is tes ted ag ainst sum rules. 



It turns out that Eqs.(4.22a) and (4.22b) are s trictly ful- 
filled by the self-consistency equations ( 4.18 ) , whereas 
the total /-spectral weight comes out slightly too large, 
J depf(e) = 3 / dep (e) — 1 . Nevertheless, the numer- 
ical difference to Eq. ( 4.22c| ) is small in the deep Kondo 



regime, where < 1/2 



V. NUMERICAL RESULTS AND RELATION TO 
THE NCA 

This Section is devoted to numerical results and their 
interpretation. In addition connections of the SCTA with 
the Non-Crossing Approximations (NCA) are discussed. 

The SCTA-Equations ( 4.18 ) have been solved numer- 
ically at the real axis z — ► u) + i0+ through iteration of 



fp (u>) . After convergence, the proper normalization of p 



and p n has been checked. Results shown here correspond 
to a flat conduction-band DOS p c {ui) = Q(D - \w\)/2D 
with cut-off D and parameters D — 10. OA , a? — —3. OA , 
Nj = 2 in the Kondo regime, where A = irV 2 p c (0) 
denotes the energy scale of the-Jipn-interacting impu- 
rity model. The Kondo energyoO takes a value Tk ~ 
0.016A. In Fig. |l0| the spectra p? and ff 1 are displayed 
for a temperature ksT — Q.ITk ■ A threshold behav- 
ior, which has been shown above to be a feature of 
these spectra at T — > (see also Appendix ^) shows 
up in the form of a sharp resonance near the bare /- 

f =f =n 

level e J . The Proj ected Spectra p , p are obtained 

by iteration of Eqs.(4.21 ) and have been normalized ac- 

cordi ng to Eq.(4.22c) ( a com mon prefactor is arbitrary in 

Eqs.( 4.21 )). From Eq.( 4.19 ) then follows the /-excitation 



spec trum ty . Since it does not strictly fulfill its sum rule 



Eq.( 4.22d ) , the latter has been used to enforce the cor- 
rect norm of p* . 

In Fig. |ll|the calculated /-spectrum is shown, together 
with an NCA result for the same set of parameters and 
temperature. It shows the well known features, a 'charge 
resonance' (CR) near the /-level and a sharp Abrikosov- 
Suhl resonance (ASR) near the Fermi edge u> = . In 
SCTA the CR is located slightly above e f = -3.0A, 
which is in contradiction to intuition (the gain in ki- 
netic energy due to hybridization should lower the /- 
level) and the NCA result. Obviously, in SCTA some 
spectral weight is transfered to energies near the Fermi 
level, and the amplitude of the ASR is in general toa 
high. As a consequence local FermiJiquid propertiesEZI 
and the Friedel-Langreth sum-ruleL3 are not fulfilled. 
The /-valence Nj(nl) as derived from Eq.(4.22d) takes a 
value « 0.66 too small for the Kondo regime. As an ori- 
gin of the misplaced spectral weight in the /-excitation 
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spectrum pf we might take the representation of the f- 
Green's function in SCTA, i.e. the self energy E cc as dis- 
played in Fig, jq (bo ttom), since it violates the normaliza- 
tion rule Eq.([L22d) . On the other hand, the character- 
istic shape of the 'ch arg e resonance' is already obvious 
in ff 1 shown in Fig. [L0|, which results solely from E-^ 
(depicted in Fig , g] t op) and enters pf through the con- 
volution in Eq. (|47l9|) . Also the /-valence is computed 



directly from p 



via the sum rule Eq. 
either. 



22d) and does 



not involve E c 

Some insight may be gained from a certain non- 
systematic limit of the SCTA, which reproduces the 
NCA : If the ladder structure in f or E cc (see Fig. §) is 
read from the left to the right, the 'ladder element' 7r (dis- 
played in Fig. H left) could be viewed as two holes of oppo- 
site spin coexisting in the /-orbital. This empty-orbital 
state then would be renormalized through intermediate 
particle-hole excited states with single /-occupancy, rep- 
resented by a (see Fig. ^ right). Following this interpre- 
tation, our bosonic Green's function II(z) appears similar 
to the NCA-propagator Gq(z) for the empty /-orbitaO . 
In fact, if the renormalization of lines within tt is ignored 
(NCA Limit), 



7r(z) -> 

the NCA equations 
of z — > z' = 2ef — z . 



-NCA 



(z) := l/[z-2ef] 



3 are reproduced after a re-definition 
Although the interpretation outlined 
above seems obvious, it may only serve as an instruc- 
tive starting point, because the Green's function G' (the 
dashed double line in diagrams) does not represent any 
physical particle. The 'formal meaning 'ofbothlland G f 
has to be viewed quite analogously to unphysical propa- 
gators Go and G\ of Resolvent-Perturbation Theory. In 
the limit of zero temperature this analogy even turns to 
an equality, as is argued in Appendix [D]. Nevertheless, 
in general both approaches are quite different, which is 
apparent if the d ef initio ns of G reen's functions and spec - 
tra in Eqs.(^2a|), ( p7|) , ( |Oc| ) and sum rules Eqs.fljLgf) 
are compared to their respective counterparts of the Re- 
solvent Method (see e.g. Ref. f79|). Most striking are the 
completely different diagram rules and the explicit use of 
the Kondo regime within the approach presented here. 

Utilizing the NCA Limit the unexpected shape of the 
'charge resonance' (CR) in SCTA is traced back to the 
self energy E-^ : In the NCA Limit the 'ladder ele- 
ment' 7? shows an imaginary part —\uit: nca {lo + i0+) = 
it5(uj — 2e* ) sharply peaked around 2e^ , which leads 
to the NCA-curve in Fig. [LI]. In going to the ladder 
approximation, this is replaced by — Im7r(w + i0+) = 
tt J de p* (e)p* (lu — e) . Accordingly — Im7r has a broad- 
ened shape like pf displayed in Fig. [h] (full line), with the 
peak located above 2e^ . Hence the CR in p' is shifted to 
higher frequencies and gets the asymmetric shape visible 
in the SCTA-curve in Fig. 0. A few more conclusions 
are drawn in the Summary. 



VI. APPLICATION TO THE LATTICE MODEL 

A. Elementary Perturbation Expansion 

The strong-coupling-perturbation expansion for the 
Anderson-lattice model, as the generic model for e.g. 
Heavy-Fermion Systems is discussed from the Feynman- 
diagram technique's point of view. The main difficulty 
encountered when performing an expansion around the 
Atomic Limit of lattice models is the so-called Excluded- 
Volume ProblernEJ (EVP), which hinders the application 
of the Linked-Cluster Theorem: Within the Resolvent- 
Perturbation formalism all parts of a lattice-diagram 
which involve a certaitu /-orbital site Ri = i are nec- 
essarily linked togetherEj . Accordingly diagrams have, to 
be broken up ('decoupled') into unlinked diagrams! 3 Hr 1 ! in 
order to obtain an expansion of thermal densities like the 
(site dependent) lattice /-Green's function F^ 1 through 
connected diagrams, leading to the known LNCAc3l3 
and XNCAE3 schemes. Specifying quasi particle inter.ac- 
tions (cumulants) consistently remains a difficult taskE3 . 
We consider here the lattice model in a fashion analogous 
to the impurity treatment in Sect. [y. It shall turn out, 
although Feynman-diagrams are obtained, that the EVP 
cannot be avoided completely; nevertheless it takes a 
much 'weaker' form than in Resolvent-Perturbation The- 
ory. 

After having decoupled states with at least one 
doubly occupied /-orbital from those with no dou- 
ble /-occupancy via a canonical transformation of the 
Anderson-lattice Hamiltonian in leading orde r V/U , we 
arrive at an expression for Fj"^ similar to Eq.(2.1C) . For 
large U , i.e. J — it reads in the interaction picture 



Fl 



x 



'.at / l\ ( rylatlO) lrylat X 

jA T > T > = ~ \ Z l Z j 

(r{cxp (-j[^£tf<n ;j (r)J 

{k- a ft-J t a)(T){f^ a f}_jl){T')) 



(6.1) 



x n ^ u <^ la!i '" 



with Z lat ^ = Tr[cxp(- /3 H lat ^)}. The unperturbed 
(atomic) lattice Hamiltonian is split into the bilinear part 



fl ,ot(0) = E^U 

k.a 



and the local /-repulsion ~ U . The latter does not influ- 
ence the time dependent operators occurring in Eq.([0]) 
and appears only in the product of e-fact ors. The per- 
turbations HQ Vi are generalizations of Eq. ( pTlc ) for each 
site i = 1, . . . , Nbz ■ Due to the fact that states with 
double local /-occupancy do contribute to the trace for 
all lattice sites p different from i and j , the e-factors in 
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general cannot be set equal to unity, in contrast to the 
single site i = j = 1 problem. The contribution of or- 
der V n to the perturbation series arising from Eq.(6.1) 
is proportional to 



(T{H^ n . ..H& ;in ■ (..z..)(r) • (..j..)(r')})' Qt(0) x 



(6.2) 



where time variables to be integrated out are omit- 
ted. Site indices i\, . . . , i n ; i, j are kept fixed for the 
moment and are assumed to cover a certain set L — 
{vi, i>2, ■ ■ ■ , v{\ of I < (n + 2) mutu al di stinct lattice sites. 
Accordingly the e-factors from Eq. ( |6.l[ ) for all sites JI € L 
act as unity since double /-occupancy is projected out 
by H^ — , and the expression takes the factorized form 



Eq. ( |6.2| ) . Wick 's theorem is now applicable to the first 
term in Eq. ( |6.2| ) , whereas the second term yields an ad- 
ditional c-number, U^l(- ■ •} i ° t( ° ) = (E ff e^) ( "" !) for 
ksT <C |£ | • That is, a bookkeeping of the number I 
of distinct sites involved has to take place when sums 
over site indices are performed. This 'weak form' of the 
Excluded- Volume Problem prevents the conventional di- 
rect cancelation of disconnected Feynman diagrams in 
the lattice /-Green's function. Nevertheless an expansion 
in Feynman diagrams is gained, i.e. any kind of 'decou- 
pling' of lattice diagrams becomes un neces sary. Possible 
applications will be indicated in Sect. VII . 



B. Local Approximation 

A from the outset more approximative approach to 
lattice models is the well known Local Approxima- 
tion (also referred to as Dynamical Mean-Field The- 
ory), which emerges from considerations of the limit 
of infinite spatial dimensiorO'H (d — ► oo) via weak- 
coupling-perturbational and -functional methods (see 
e.g. ReffJ§4jj93) as well as strong-coupling-perturbation 
theoriesEHHH^ doing with Hubbard cumulantsEj . For 
the AnderjsapJattice model the different— d-. — > oo- 
approachesOo'Ell lead to identical resultsOEJ , which 
becomes obvious if self-energies are eliminated from 
equations in Refs. |S7| , |93| . The lattice /-Green's func- 
tion F 1 ^ in. fc-space turns out equivalent to the XNCA 
expression^ , 



*kcr 



(F a (z))-i-V* 



1 



with an effective 

impurity (EI) model's conduction-electron Green's func- 
tion Q c a determined through the lattice-self-consistency 
condition (Lattice-SC) , ^^^'(z) = F a {z) . The 
EI model /-Green's function is designated as F a . In cal- 
culating F a numerically via the SCTA theory within a 



loop that iterates Q c a to convergence, the approach de- 
scribed in previous Sections is extended to the U — > oo 
Anderson-lattice model. Unfortunately, we encounter a 
difficulty well known from the XNCAE3 : Except at ex- 
tremely high temperatures the spectrum of Q% develops 
a peak structure near the Fermi level, indicating an en- 
hancement of the Kondo effectCJ ; but if Q c a is re-inserted 
into the SCTA , the spectrum of F a does not alter sig- 
nificantly, i.e. it does not show any tendency toward a 
hybridization gap in the Abrikosov-Suhl resonance. Con- 
sequently the Lattice-SC cannot be fulfilled. Comments 
on this result are given in the following Summary. 



VII. SUMMARY AND CONCLUSIONS 

A detailed description of a self-consistent strong- 
coupling-perturbation theory for the Anderson model has 
been given, which jexlonds a Feynman-diagram technique 
proposed recentljH't3 . It was shown in Sect, [n] that 
Wick's theorem can be applied directly within a perturba- 
tion expansion in the hybridization V , i.e. the coupling of 
a localized orbital (/-orbital) to delocalized conduction- 
electron states. Although the strong Coulomb inter- 
action U of electrons on the /-orbital is included in 
the unperturbed Atomic Limit, conventional Feynman- 
diagram rules are obtained in the sequel, involving the 
two-particle vertices ~ V and ~ V 2 /U shown in Fig. |l]. 
This kind of Feynman-perturbation theory starts from a 
canonical transformation of the Hamiltonian, which elim- 
inates in leading order V/U the charge fluctuations into 
states with double /-occupancy. It has also been pointed 
out in Sect. [lj that an application to the (multi channel) 
Kondo model is apparent. 

Based on the diagram rules for U — > oo the Self- 
Consistent T-Approximation (SCTA) has been derived 
in Sect. 



Ill . By means of irreducible vertex functions 



compliant to symmetry relations (Crossing Relations) a 
matrix-type self energy S has been set up; it is displayed 
in Figs. || , [6|. This skeleton-self energy constitutes a <&- 
derivable approximation for a set of unphysical propaga- 
tors as well as a prescription for constructing the physical 
/-Green's function from these propagators. It contains 
ladder sums of the diagram elements shown in Fig. ^| to 
infinite order V ; the latter have been identified as impor- 
tant for the formation of the Kondo effect. A version of 
the SCTA suitable for the Kondo regime emerged in Sect. 
IV from an expansion of S in a non-perturbative small 
parameter. Only one of the unphysical Green's functions 
remains renormalized in this case, r eprese nted through 
the self-consistent SCTA-Equations (4.18). Some rela- 
tions to the Non-Crossing Approximation (NCA) were 
explored in Sect. Furthermore the SCTA contains 
the NCA as a certain non-systematic limit. 

The numerical results for the physical /-Green's 
function presented in Sect. reproduce the known 
correlation-induced features of the /-excitation spec- 
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trum. Concerning Fermi-liquid properties, it became ap- 
parent that the SCTA does not furnish results which im- 
prove those of the NCA . Nevertheless, the formulation of 
enhanced approximations will be straight forward since 
standard techniques for Feynman diagrams (e.g. Parquet 
equations) can be used, and the expansion for the Kondo 
regime is not restricted to specific diagram classes (see 
Appendix [bJ). Going beyond the summation of ladder 
diagrams in the ^-functional will involve 'true' retarded 
vertex functions. The importance of vep tox . co rrections to 
the NCA has already been pointed outilSESEH . Here ver- 
tex corrections to the matrix-type self energy S yield a 
consistent modification of the self-consistency equations 
as well as the representation of the physical /-Green's 
function. Actually both have significant influence on the 
/-spectrum and require corrections, as has been argued 
in Sect. fv|. Further studies based on the technique de- 
veloped here might therefore be fruitful with regard to 
dynamical properties of the spin-degenerate model in the 
problematic region of very low temperature. 

Two possible ways of extending the technique to the 
(U — > oo) Anderson-lattice model, an elementary per- 
turbation expansion on the lattice and the application 
within th e Lo cal Approx imation, have been discussed 
in Sects. VIA and VI B respectively. For the elemen- 



tary expansion it turned out that any contribution to the 
/-Green's function is indeed decomposed into Feynman 
diagrams by virtue of Wick's theorem, but in addition 
the number of mutual different /-orbital sites actually 
involved enters through a c-number ('weak' Excluded- 
Volume Problem). Thus an expansion in connected dia- 
grams may not obtained in a simple fashion. Neverthe- 
less, the definition of quasi particles and their interac- 
tions seems manageable. Accordingly the approach looks 
quite suitable for models with a small number of impuri- 
ties with local interaction, e.g. the two-impurity multi- 
channel Kondo modeS . Concerning the Anderson- 
lattice and i-J-model we believe that our method will 
also show some impact, especially on the study of non- 
local spin fluctuations via strong-coupling theory for e.g. 



High-T c -Superconductor 



jfermion Sys- 



tems near a quantum-phase transitioned 

In the Local Approximation (or D ynam ical Mean-Field 
Theory, reviewed briefly in Sect. VI B) the /-Green's 



function of the lattice model is given through that of 
an effective impurity (EI) model subject to a lattice- 
self-consistency condition (Lattice-SC). The SCTA equa- 
tions from Sect. |v| were used to solve the EI model, but 
convergence to the Lattice-SC has not been achieved; 
a feature already known from the— JNCA within Lo- 
cal Approximation (i.e. the XNCA)o. Such a behav- 
ior is not found if the Lattice-SC for the Anderson- 
lattice model is combined with methods different 
strong-coupling perturbation theory (e.g. QMC 
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i a prm alizationEHU or weak-coupling perturba- 



numerical cL-,, 

tion theoryE£3'E£El) , or if strong-coupling theory is used 
within a different Lattiee-SC for the same model, as is 
the case of the LNCAE3'E3. In order to clarify this in- 



consistent situation we propose considering alternative 
strong-coupling expansions for the lattice model, which 
do not rely on the Local Approximation or the Resol- 
vent Method. The elementary expansion outlined in Sect. 
VIA may serve as a starting point. 
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APPENDIX A: EFFECTIVE HAMILTONIAN 



An effective Hamiltonian H" , where the mixing H\ 2 
involving local double occup ation is removed to first or- 
der V/U is derived from .Eq. (2 1 ) in two steps: The first 
unitary transformationHE23 is H — ► H' , 



H' = e iS He~ iS = H+ [iS, H] + %[iS, [iS, H}] 
with a generator iS subject to the requirement 
[iS, (H c + H loc )} = -HY 2 . 



An appropriate choice is 



iS = -= — V* J a {k) n f _ a flc k0 



h.c. 



(Al) 



with J a (k) = V k /(el+U-£ k a) ■ Since is kept in H' , 
additional operators mixing Hubert space sectors 7io,i 
and H.2 emerge through [iS, H^] , which is of order V/U . 
To this order they vanish in a second transformation 



H' -> H" =H'+ [iS',H'\ + . 
using a generator iS' which fulfills 

[iS',(H c + H l ° c )] = -[iS,HV] 
It is explicitly given by 



Cfco 



h.c. 



BZ 



k.q.cr 



involving coefficients K. a (k,q) of order V 2 /U 2 . The 
resulting effective Hamiltonian H" reads to first order 
V/U, 



H" = H c + H l 
where the Coulomb repulsion U in H loc is slightly shifted, 
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U = U +T7—Y 



N B z ^ e J a + U - e k<J 



and a spin-spin interaction emerges, which is part of 



N BZ ^ 2 



(/<r<4 -vf-vCkv + ftcl,-a C k-<rfo) + h.t 



The exchange integral originating from Eq.(Al) is sim- 
plified for large U and V k = V to 

JAk)V q *^ \V\ 2 /(e f + U) = - J v , 



and H' takes the form Eq.(2.4) . 



APPENDIX B: THE LOOP THEOREM 



In Section IV counting of a certain kind of loops is 



used to separate out important self-energy diagrams for 
the Kondo regime. The underlying 'loop theorem' for 
skeleton diagrams is derived as follows. Consider a closed 
path in an arbitrary diagram, as introduced in Section 
IV below Eq. fl4.7| ) : It consists solely of equally directed 
renormalized Green's functions G* (dashed double lines, 
see Fig. |^) and vertices as displayed in Fig. [I] (a) . Vertices 
may be passed from any incoming dashed line to any 
outgoing dashed line, i.e. the path may involve both spin 
directions. None of the diagram's lines and vertices is to 
be touched more than once. In Fig. ^ (a) some examples 
are shown. The analytical contribution from a closed path 
running on k dashed double lines involves also k distinct 
vertices and is proportional to 



1 



(Bl) 



±- 



G f ( 



iv x ) 



,G f (i 



It is coupled to the remainder of the diagram through k 
fixed bosonic frequencies labeled V\ , . . . , v k . These are 
kept unequal in-pairs through suitably fixed infinitesi- 
mal incremental^ v\ — > v\ + Si . The sum over the free 
fermionic loop- frequency uj n is rewritten as usualo , 



L {k) =T J de 1 ...e k p f (e 1 )...p f (e k ) 

[!-/(*)] 



dz 



2ni (z - ivi - e\) . . . (z - iv k - e k ) 



where is expressed via its spectrum Eq. ( |4.6| ) . The 
contour Y encircles all zeros of the denominator. The 
term [1 — Hz)} is a consequence of the exponential factor 
in Eq.(Bl), which keeps track of non-normal ordering 
of operators (ref. Sect. II B). Performing the contour 
integral leads to 



L (fc) 



Now all other free frequencies in the diagram are in- 
tegrated out one after the other, and the diagram ap- 
pears as a sum over k terms. These are regular for fi- 
nite temperature, and the infinitesimal increments intro- 
duced above become unnecessary. The function p* (e) [1 — 
/(e)] ~ k discussed in Section IV shows up in each term 



here and gives the diagram the small overall order of mag- 
nitude k — exp(/3(e^ + Tk)) in the Kondo regime. The 
argument becomes apparently invalid if at least one of the 
Green's functions & on the path is reversed in orienta- 



tion, i.e. 



-e) , or is replaced by a conduction- 



electron propagator pf(e) — > p c (e) ■ In a straight-forward 
generalization, the resulting 'loop theorem' is stated: p 
non-overlapping closed directed paths in a skeleton dia- 
gram are sufficient for an order of magnitude ~ (k) p . 

As an application we determine the order of magni- 
tude of self energies T, aa ,a'a = ff,cf,fc,cc, without 
reference to any specific approximation. The result has 
been quoted in Eq. ( [Qj} ) . Consider at first a diagram to 
S cc , or more generally an arbitrary connected dressed 
skeleton diagram with one incoming and one outgoing 
conduction-electron line G c . It shows necessarily at least 
one closed path as defined above: We construct a path 
running solely on dashed lines & , starting on an outgo- 
ing dashed line of an arbitrary chosen vertex within the 
diagram. We always follow the direction of arrows, and 
no line or vertex is visited twice. Owing to the fact that 
a vertex (displayed in Fig. |l](a) left or right) entered on 
the way can be left via at least one dashed line, the path 
cannot terminate within the diagram. Therefore (i) the 
procedure can lead back to the vertex we started from, 
and a closed path is completed; (ii) a vertex with two 
incoming dashed lines (Fig. [l] (a) right) is re-entered for 
the first time on its hitherto unused line, which consti- 
tutes a closed path, too; (iii) the path we are constructing 
leaves the diagram via an external outgoing dashed line, 
if present. Possibility (iii) is ruled out for the diagrams 
to E cc , therefore these contain at least one closed path 
which contributes a factor k to the self energy's order of 
magnitude. Diagrams showing exactly one such path are 
easily found (see Fig. || bottom) , and it follows that in 
general E cc ~ (k) 1 . In quite the same way it is seen that 
~ (k) 1 and ~ (k) 1 . For the latter a closed path 
is constructed starting on an incoming line of an arbi- 
trary chosen internal vertex and proceeding in opposite 
direction to arrows. The situation is altered if T,^ is con- 
sidered: Here a path can leave the diagram (case (iii)), 
that is may contain diagrams of order ~ (k)° = 1 . 
Actually such diagrams do exist (e.g. those shown in Fig. 
|top), thus E^ - (k)°. 
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APPENDIX C: EQUATIONS FOR PROJECTED 
SPECTRA 



A set of self-consistency equations is derived for the 
Projected Spectra introduced in Section IV . According 



to the definition of p in Eq.( 4.20a) we consider 
pf(w)[l-f(u)] = 

&{u)) Vi)lm£''V + zO + )[l-/H] , 



where Eq.(4.18b) has been used with z = uj + i + . T he 
self energy's imaginary part is taken from Eq.( 4.15| ) at 

iuji — ► uj + z0+ , 



-/(«)] = Gf(oj) 



(CI) 



xV 2 / dep c (e)p n (e + uj)[~f(e)~g(e + uj)]f(-uj) 



with 1 — f(uj) = f(—u>) . For further transformation of 
expressions like this, a couple of relations involving Bose 
and Fermi functions (designated as g and /) is useful, 



g(x-y)=f(x)[l-f(y)}/[f(y)-f(x)} , 
f(x-y)=f(x)[l+g(y)}/[g(y) + f(x)} . 



(C2a) 
(C2b) 



For use in Eq.(Cl) the former is applied to g(e + uj) in 
the first row of the right hand side in 



.]/(-") 



[-/(e) -g{e + w)][f{-w)- /( £ ) 

+ [-/(e)-. 9 (e + c)]/(e) 
/(e)[-l-.g(e + ^)] . 



Thereby Eq.(Cl) becomes 



(?(u)[l-f(w)]= G / (w) 



xF 2 / dep c (e)/(e)p n (e + ^)hl-. 9 (e + W ) 



(C3) 



By multiplying Z/Z to both sides Equation ( 4.21a ) is 
obtained immediately. 



Its counterpart Eq. (4.21b) is derived in a similar fash- 
ion: With Eq.(4.11) the spectrum p u (ui) is written 



P n = -iim n = 



n 



2 ( — Im (tt)/tt 



Ilm (cF) 



7r and cr are taken directly from the graphical definition 
shown in Fig. , 



7r(z) 



(C4a) 



a{z) = NjV 2 ( (dede 1 p c {e)p {e 1 ) 1 -^ — £i-£i . (C4b) 



z — e — e' 



The e ffect of a prefactor [—1 — g(w)] to Eqs.( C4a ) and 
(|C4b ) is made explicit by application of Eqs.( |C2[ ) Col- 
lecting the resulting parts together, we get 



j5»[-l- 5 (o0 



n( w ) 



x ds 



P{e)[l-f{e)] 



+ NjV 2 p c (e)[l-f(s)} 



{ l*Mr 

xpf(oj-e)[l-f(w-e)] 



Since | tt | 2 3> K the first term under the braces may be ne- 
glected with Eq.(4.7) , and Eq.(4.21b) appears in writing 



Z/Z on both sides of the above formula. Furt herm ore, 
the result confirms the order of magnitude Eq.( |4.16| ) di- 
rectly for the special case of the ladder approximation. 



APPENDIX D: SOME EXACT RELATIONS 



1. Sum Rules 



Exact sum rules are proven for spectra p 



=f 



=n . f 

p , and the physical /-spectrum p 1 . These are ver- 
ified within the Self-Consistent T- Approximation. Re- 
sults have already been stated in Section IV through 
Eqs.(4.22). The spin cr = ±1 is arbitrary in the fol- 
lowing, since no external magnetic field is considered. 

The total sp ectra l weight of the propagator is de- 
rived from Eq.( 3.2a ) as ({f a , /j}) = 1 , accordingly the 
correspondin g spec trum defined in Eq.(^^) fulfills the 
sum rule Eq. ( 4.22a ) . In the same fashion, the spectral 



weight of II follows with Eqs.(4.17), (4.14) in the form 



J de p n {e) = 6(0_) - n(0+) 
= (»{«{>- <(1 - uf)(l - n{))= 5>#- 1 



and by intro ducing jo* via (n£) = J dep> (s)f(e) , the 
sum rule Eq.(4.22b) is confirmed. 

In order to derive a sum rule for the Projected Spectra 
we first note the physical /-spectral-weight factor per 
spin, 



(l-nl)=-(Z/Z)Gl(0 
Z 



(Dl) 



dspf(e)[l-f(e))= I de-p\e) 



Consider also the probability of the local /-orbital being 
empty, 



<(l-nf)(l-n{)) = (Z/Z)II(0 + ) 



f\\ - 



(D2) 



| f dep n (s)[-l-g(e)}= / def(e) 
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In the last two expressions the respective first equality 
sign stems from the fact that states with doubly occupied 
/-orbital do not con tribu te to the traces (this has been 
explained in Section 



II B 



). At U — > oo the probability 
of double /-occupancy is zero, and with the probabil- 
ity (n£(l — n>L a )) for single /-occupancy a completeness 
relations holds, 

((l-n{)(l-n{))+5>£(l-»^)> = l . 



filled Fermi sea of conduction electrons. It is unaf- 



fected by the perturbation in Eq.( 2.11 ) , and its energy 
is significantly lower than the ground 
state energy Eg ~ E c + — Tk of the 'physical' sub- 



space Ho,i with /-occupancy < 1 ; see also Sect. IV A 
Since /]|2;FS C ) = and / CT |2;FS C ) = ±| - cr; FS C ) 



follows pl^> (e) = and 



[%\-a;FS c )\U(e-[2ef +E c -E i }) 



This lead s with Eqs. (|Dl| ) and (|D2D to the sum rule 
Eq.( f4.22c| ) stated in Section |Tv|. 

The physica l /-ex citation spectrum is normalized ac- 
cording to Eq.( 2.10| ) , i.e. 



J dep f (e) = F (T (0^)-F rr (0 + ) = (1-nt, 



By use of Eq. (|Dl| ) , the sum rule Eq.( 4.22c| ) is evident. 

A verification of sum rules in SCTA is performed for 
pf and p n via the behavior of corresponding propagators 
at infinity, 



1*1 



oo 



G f (z) -> a f /z , n(z) -> a u /z 



From the self-consistency equations (4.18) at \z\ 
get 



oo we 



1 , a 1 



dspf(e)[2f(e)-l] 



i .e. th e values from exact sum rules Eqs. ( 4.22a ) and 
( 4.22b] ) are reproduced. The /-spectrum's norm alization 
in SCTA is obtained via integration of Eq.( |4.19| ) , 



dw p (lo) 



dhj p\uj) J dep n (e) + j dujp n {uj) 
:>, I dw? f (u) - 1 . 



Corrections ~ (k) 1 from application of the sum rule 



Eq.(4.22b) are negligible in the Kondo regime. O bvious ly 
the norm of p' does not match the sum rule Eq. ( 4.22d ) . 



2. Spectral Decomposition at Zero Temperature 

Consider the unphysical spectra ~p> and p n at zero tem- 
perature T = . A spectral decomposition of G£ leads 
to pl(e) = pl i+) (e) + p f J~\e) where 



c*)=El< 



S(eT[E a -E }) 



at T — > . An /-creation operator f\ is used i n p 



•/(+) 



The ground state |\l/o) of the Hamiltonian Eq.(2.11 



the full Fock space of /- and conduction electrons is 
|*o) = |2;FS C ), a doubly occupied /-orbital plus the 



From the complete set {I^q)} of eigenstates only the 
subspace {l^i)} = Ho.i contributes to the overlap in- 
tegral, with energies Ei > Eq . Accordingly the spec- 
trum pl(e) ~ Q(E Th — e) shows a threshold at energy 
E Th = 2£ f — Eq + E c slightly above the bare /-level 
. The same threshold emerges in /^(e) ~ Q(E Th — e) : 
From a similar analysis p^{e) is obtained in the form of 
p f with | - cr; FS C ) replaced by |0;FS c ) . At zero tem- 
perature the unphysical spectra p£ and p 11 match the so- 
called 'jepc spectra' p a and po of Resolvent-Perturbation 
TheoryEj via 



Pa (e) = pl{2ef - e) , Po (e) = ^(2^ - e) 



(D3) 



Accordingly an algebraic singularity is expected in p^ n 
at energy approaching E Th from below. It ha s ne verthe- 
less to be noted that the strict equality Eq.(|D3|) is lost 
at any finite temperature and is likely to be lost if ap- 
proximations are cons ider ed (e.g. SCTA vs. NCA). Also 
is a relation like Eq.(D3) not found for the Projected 
Spectraontroduced in Eq.( 4.20| ) and the so-called Defect 
Spectrally (also called e.g. 'negative-frequency spectra'Ea) 
of NCA theory. 
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FIG. 1. Vertices of the Feynman-diagram expansion. 
Das hed lines correspond to local Green's functions Gr '- ' , 
Eq.( 2.16 ), full line s stand for conduction-electron propaga- 
tors G c , Eq. (|2 . 1 7|) . Vertices (a) originate from (see 



Eq.(2.1c)) and are interpreted in the subspace without doubly 
occupied local /-orbital: A transfer band- to /-electron (left) 
or vice versa (right) is permitted only if a /-hole of opposite 
spin is present. In ( b) t he spin- flip and charge scattering 
vertex of H J (see Eq.(2.4)) are shown on the left and right, 
respectively. 




FIG. 2. Simple ladder diagrams from particle-particle 
(pp) and particle— hole (ph,ph) vertex contributions to the 
/-Green's function at U — ► oo . These may show logarith- 
mic divergences, see text. 
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FIG. 3. (a): Notation of renormalized Green's functions 
occurring in skeleton diagrams of the U — * oo problem. A def- 
inition is given in Eq. ( ft.2| ) . (b) : 2 x 2-matrix-Green's function 
(left), its matrix elements are shown in (a) , and correspond- 
ing bare propagator (ri ght ) . (c): Dyson's equation involving 
matrix-self energy Eq.(ti.4|) . 
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FIG. 6. Self energy in Self-Consistent T- Approximation. 
With spin restriction at vertices, all exchange parts present 
in the vertex function in Fig. ^ do not contribute to S and 
internal spin sums do not occur. The self energy is exact to 
V 2 (first row), higher orders separate into ladder sums from 
(pp)-channel (2nd row) and (ph)- (3rd row) and (pfr)-channels 
(4th row). 



FIG. 4. (a): Vertex used with matrix-Green's func- 
tion from Fig. H (b) for U — > oo . It collects both hy- 
bridization types shown in Fig. |l| (a) . 'Spinor' indices 
{ j() ) J,C ) = lO^' indicate canonical operators l' ' = f 
2' ' = . (b): Vertex combination used in expansion of 
matrix-self energy £ . 
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FIG. 5. Exact skeleton representation of matrix-self en- 
ergy at U — > 00, with full Green's function and vertex from 
Figs. ^ (b) and ^ (b) . Due to the vertex' spin restriction no 
'Fock diagram' is present. All higher orders are collectively 
expressed through the irreducible vertex function Y 2 P . 




(a) 



(b) 



FIG. 7. Examples of diagrams to which do not con- 
tribute in the Kondo regime: (a): due to the loop theorem 
(possible closed paths are indicated by dotted lines); (b): due 
to the presence of G fc or G cf . Notation follows Fig. | (a) . 
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FIG. 8. Self energy S in the Kondo regime, consisting 
only of (pp)-type diagrams with Green's functions Cr (dashed 
double line) and bare conduction electrons G c (full lines). 
Off-diagonal elements E c ^ and E^ c as well as all (ph)- and 
(p/i)-type diagrams included in Fig. ^ are negligible in the 
Kondo regime, examples are shown in Fig. [?]. In each pla- 
quette indicated by a small double arrow, parallel band- and 
/-lines may be interchanged, leading to the full set of dia- 
grams to be considered. 
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FIG. 9. Definition of 'ladder elements' occurring in Fig 
. Spin degeneracy leads to the prefactor Nj = 2 in a(ivk) ■ 
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FIG. 11. Local /-excitation spectrum of the U — > oo 
Anderson model from Self-Consistent T-Approximation (full 
line) and NCA (dashed dotted line). Inset: Abrikosov-Suhl 
resonance near the Fermi energy uj = . Parameters as in 
Fig. 0. 
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FIG. 10. Spectra (full line) and /5 11 (dashed line) from 
numerical solution of SCTA equations (see text). Parameters 
are e ! = -3.0A, D = 10. OA , Nj = 2 , and k B T = 0.12k . 
The spectral weight turns to zero for energies above the 
threshold. Inset: Neighborhood of the threshold (note the 
different vertical scale). 
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